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ABSTRACT 

We present C/iandra X-ray imaging and spectroscopy for 14 quasars in spatially resolved pairs tar- 
geted as part of a complete sample of binary quasars with small transverse separations drawn from 
Sloan Digital Sky Survey (DR6) photometry. We measure the X-ray properties of all 14 QSOs, and 
study the distribution of X-ray and optical-to-X-ray power-law indices in these binary quasars. We 
find no significant difference when compared with large control samples of isolated quasars, true even 
for SDSS J1254+0846, discussed in detail in a companion paper, which clearly inhabits an ongoing, 
pre-coalescence galaxy merger showing obvious tidal tails. We present infrared photometry from our 
observations with SWIRC at the MMT, and from the WISE Preliminary Data Release, and fit simple 
spectral energy distributions to all 14 QSOs. We find preliminary evidence that substantial contri- 
butions from star formation are required, but possibly no more so than for isolated X-ray-detected 
QSOs. Sensitive searches of the X-ray images for extended emission, and the optical images for optical 
galaxy excess show that these binary QSOs — expected to occur in strong peaks of the dark matter 
distribution — are not preferentially found in rich cluster environments. While larger binary QSO 
samples with richer far-IR and sub-millimeter multiwavelength data might better reveal signatures of 
merging and triggering, optical color-selection of QSO pairs may be biased against such signatures. 
X-ray and/or variability selection of QSO pairs, while challenging, should be attempted. We present 
in our Appendix a primer on X-ray fiux and luminosity calculations. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei - 
quasars: emission lines 
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1. INTRODUCTION 

Supermassive black holes (SMBH; > Kf' Mq) mainly 
grow from gas accretion (Soltan 1982; Merloni & Heinz 
2008) in the cores of active galaxies. Luminous quasars, 
which host SMBH up to about a billion solar masses, are 
already in place by redshift z^7 when the universe was 
much less than a billion years old (Fan et al. 2001, 2006; 
Mortlock et al. 2011). While for decades SMBH accretion 
was thought of largely as a product of galaxy evolution, 
it is also seen now as a principal driver of that evolution, 
based in large part on the surprisingly tight correlation 
between the mass of SMBH and the mass or velocity dis- 
persion of their host galaxy bulges (Ferrarese & Merritt 
2000; Giiltekin et al. 2009). If SMBH are the product 
of a sequence of galaxy merger episodes (e.g., Hernquist 
1989; Kauffmann & Haehnelt 2000; Hopkins et al. 2008), 
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then binary SMBH are an inevitable outcome of galaxy 
assembly. As we seek greater understanding of the cos- 
mic co-evolution of galaxies and SMBH, the relatively 
rare binary SMBH hold unique interest and promise, and 
particularly so for the most luminous examples. 

Luminous quasars have always inhabited a relatively 
small fraction of galaxies. Studies of the clustering prop- 
erties of quasars (two-point correlation functions) indi- 
cate that the bias of quasars relative to underlying dark 
matter increases rapidly with redshift, implying that 
quasars inhabit rare massive dark matter halos of sim- 
ilar mass (Mhaio ^ 10^^ h'^Mo) at every cosmic epoch 
(Porciani et al. 2004; Groom et al. 2005; Myers et al. 
2006, 2007a; Shen et al. 2007; Ross et al. 2009). Though 
quasars with a close (< 1 Mpc) quasar companion at 
comparable optical luminosity constitute only ~ 0.1% 
of optically selected quasars overall, that represents a 
strong excess at small separations over the extrapolation 
from the larger scale QSO spatial correlation function 
(Hennawi et al. 2006; Myers et al. 2008; Hennawi et al. 
2010; Shen et al. 2010a). Indeed, the surprisingly large 
number of binary quasars in the universe (Djorgovski et 
al. 1987; Myers et al. 2007b; Hennawi et al. 2006) is a key 
underpinning of the merger hypothesis. Multiple authors 
(Djorgovski 1991; Kochanek et al. 1999; Mortlock et al. 
1999; Myers et al. 2007b) have noted that an excess of 
binary quasars could be due to tidal forces in dissipative 
mergers that trigger inflow of gas towards the nuclear 
region — and hence strong accretion activity in the nu- 
clei of merging galaxies. However, this picture remains 
incomplete; are mergers the cause of the observed ex- 
cess of binary quasars or rather is the excess of binary 
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quasars the result of enhanced smah-scale clustering for 
the merger-prone halos that host quasars? The measured 
small-scale excess including binary quasars {R > 100 kpc) 
may not be due to mutual triggering, but rather simply a 
statistically predictable consequence of overdense, group- 
scale environments (Hopkins et al. 2008). These contro- 
versies motivate detailed studies of binary quasars — and 
binary Active Galactic Nuclei (AGN) in general. 

At high redshifts, where the merging process is likely to 
be efficient (e.g., Springel 2005), binary AGN are difficult 
to resolve. At more recent epochs, where they could be 
resolved, the merger rate is lower (Hopkins et al. 2008). 
Nearby examples exist, however. The merger hypothe- 
sis is supported by (1) the existence of spatially-resolved 
binary active galactic nuclei in a few z < 0.1 galaxies 
with one or both of the nuclei heavily obscured in X- 
rays (NGC6240; Komossa et al. 2003; Arp299, Zezas et 
al. 2003; Mrk463, Bianchi et al. 2008), by (2) the un- 
usual BL Lac-type object OJ 287 (Sillanpaa et al. 1988; 
Valtonen et al. 2011) whose binary nature is still under 
considerable debate (Villforth et al. 2010), and perhaps 
by (3) X-shaped morphology in radio galaxies (e.g., Mer- 
ritt & Ekers 2002; Liu 2004; Cheung 2007). In addition, 
CXOC J100043.1+020637 contains two AGN resolved at 
0.5" (~2.5kpc) separation in HST/ACS imaging, which 
have a radial velocity difference of A?; = 150 km/s, and 
appear to be hosted by a galaxy with a tidal tail (Com- 
erford et al. 2009; Civano et al.' 2010). 

A recent flurry of searches for candidate close binary 
AGN (with sub-kpc projected separations) has mostly 
involved spectroscopic (unresolved) binaries. Some show 
both broad and narrow emission lines with significant 
velocity offsets, such as SDSS J153636. 22+1044127.0 
(Boroson & Lauer 2009) or SDSS J105041.35+345631.3 
(Shields et al. 2009). Some may be true binary SMBH. 
SDSS J092712.65-F294344.0 may be a binary SMBH with 
a single disk (Bogdanovic et al. 2009; Dotti et al. 

2009) , or a single SMBH that has been "kicked" due to 
anisotropic emission of gravitational radiation near coa- 
lescence (Komossa et al. 2008) . Some may be similar to 
spatially unresolved quasars with double-peaked broad 
emission lines (e.g., Strateva et al. 2003). Debate sur- 
rounding the various interpretations persists (e.g., Lauer 
& Boroson 2009; Wrobel & Laor 2009; Tang & Grindlay 
2009; Vivek et al. 2009; Chornock et al. 2010). Many 
spectroscopic binary AGN candidates with narrow emis- 
sion lines only have been selected from the Sloan Digital 
Sky Survey (SDSS; York et al. 2000) based on double- 
peaked [O III] AA4959, 5007 emission lines in their fiber 
spectra (Wang et al. 2009; Smith et al. 2010; Liu et al. 

2010) . Some remarkable objects have been found (e.g., 
Xu & Komossa 2009), but several scenarios can pro- 
duce double-peaked narrow emission lines, including pro- 
jection effects, outflows, jet-cloud interactions, special 
narrow-line region (NLR) geometries, or even a merger 
where one AGN illuminates two NLRs. Near-infrared 
(near-IR) imaging and optical slit spectroscopy can re- 
veal genuine double-nuclei (Liu et al. 2010), which con- 
stitute only about 10% of the candidates (Shen et al. 
2010b). Chandra imaging is underway now to confirm 
these nuclei as AGN by resolving their luminous hard 
X-ray emission. Fu et al. (2010) imaged 50 double- 
peaked [OIII]A5007 AGN from the SDSS with Keck-II 
laser guide star adaptive optics, confirming that most 



(70%) are probably single AGN. 

Spatially resolved, confirmed mergers of broad line 
AGN may well be the most useful systems for tracing 
the physics of the early merger process because they 
probe ongoing mergers, and because the spatial and ve- 
locity information, especially when combined with well- 
resolved spectra providing separate black hole mass es- 
timates, offer more constraints on the properties of the 
merging components and the evolution of the merger. 
Examples of resolved binary broad line AGN in con- 
firmed mergers are virtually unknown. Probably the 
best example to date is SDSS J1254+0846 (Green et al. 
2010), which clearly shows tidal tails from the ongoing 
merger. Spatially-resolved active binary mergers such as 
these provide by far the strongest constraints on merger 
physics at kiloparsec scales. Even when such obvious 
merger signatures arc not available, other probes of the 
properties of binary quasars such as their environments, 
spectral energy distributions (SEDs), and nuclear and 
host galaxy properties provide useful information to help 
distinguish which systems may be undergoing merging or 
triggering, and to elucidate merger physics itself. 

In this paper, we probe the multiwavelength proper- 
ties of a small but uniform sample of binary quasars, 
described in §2. Using Chandra X-ray imaging, in §3 
we study the high-energy SEDs of binary quasars (§ 3.1) 
and how they compare (in §3.2) to a subset of QSOs 
imaged in X-rays by Chandra as part of the Chandra 
Multiwavelength Project (Green et al. 2004; Green et al. 
2009). In §4 we present deep IR imaging we obtained 
at Mt Hopkins using the SAO Wide-field InfraRed Cam- 
era (SWIRC) on the 6.5m MMT, to further examine the 
SEDs of binary QSOs. Template-fitting to our multi- 
wavelength SEDs presented in § 5 will allow us to con- 
trast binary QSOs directly with a large sample of iso- 
lated QSOs from the Chandra Multiwavelength Project 
(ChaMP) in an upcoming paper. 

Binary quasars are expected to frequent massive dark 
matter halos, which we test in § 6. We look for evidence 
of any local hot intra-cluster medium (ICM) indicating 
a host group or cluster in §6.1. X-ray cluster detec- 
tion avoids some of the pitfalls of optical/IR selection — 
namely, projection effects and red-sequence bias towards 
evolved galaxy populations. Our X-ray imaging is sen- 
sitive even to poor clusters and groups with high M/L 
(Barkhouse et al. 2006), despite the presence of bright 
quasar point sources (Green 2005). Analysis of optical 
images of these fields is described in § 6.2. We obtained 
N0A0/4m-M0SAIC images on Kitt Peak to deeply im- 
age the quasars (§6.2.1) in a search for signs of merger 
activity or optical galaxy overdensities associated with 
a host cluster, but SDSS imaging (§ 6.2.2) largely turns 
out to be more useful than what we could obtain in vari- 
able weather. We conclude with a brief discussion of our 
findings in § 7. 

Throughout, we assume a cosmology with Cl^ = 
0.3, Oa = 0.7, and Hq = 72 km s'^ Mpc'^ 

2. BINARY QUASAR SAMPLE 

Ongoing mergers hosting two luminous AGN are rare, 
so while a handful of serendipitous examples exist, huge 
volumes of sky must be searched to find them systemat- 
ically. The SDSS provides a large sample of multicolor 
imaging and spectroscopy for this purpose. 
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Candidates in our sample (see Table 1) were drawn 
from the photometrically classified Type f (broad line) 
quasar catalog of Richards et al. (2009). All pairs of ob- 
jects with component separations of 2.9" to 6", "UVX" 
or "low-z" classification flags ^ set and g < 20.85 were 
targeted for follow-up spectroscopy. Selecting "low-z" 
and "UVX" objects produces a sample that populates 
the redshift range 0.4<z<2.4. At angular separa- 
tions of < 2.9" we supplemented our sample with SDSS 
J0740-h2926 from Hennawi et al. (2006). 

Our relatively uniform parent sample allows us to place 
these systems in their larger cosmological context, which 
is crucial if we are to understand the role of merger- 
triggered supermassive black hole accretion, and its rela- 
tionship to galaxy evolution. By selection, the two com- 
ponents of these quasar pairs are likely to have a wide 
projection on the sky, which makes them useful for pro- 
viding morphological constraints on merger models. 

As outlined in Table 2, the quasar pairs we consider 
in this work are confirmed, through spectroscopy, to be 
broad line quasar pairs with components that are prox- 
imate in velocity space ("binary quasars"). Five of our 
pairs had one object already confirmed spectroscopically 
by SDSS. (The ^^55" minimum SDSS fiber separation 
usually precludes SDSS spectra for both members.) The 
majority of the confirmed binary quasars we study were 
followed up and spectroscopically identified in previous 
works (Hennawi et al. 2006; Myers et al. 2008). Ex- 
ceptions are SDSS J0813-F5416 and SDSS J1254-F0846. 
These were spectroscopically confirmed on February 9- 
11, 2008 at Kitt Peak National Observatory using the 
R-C spectrograph on the Mayall 4-m telescope and the 
KPC-lOA grating. A 1.5" wide long-slit was oriented at 
a position angle to simultaneously observe both candi- 
dates, and 20 minutes of exposure time was sufficient to 
identify the candidates as broad line quasars (the faintest 
quasar in these two binaries is at 17 = 20.3). Spectra were 
reduced using the LowRedux^ package, better quality 
spectra of the QSOs in SDSS J1254+0846 were obtained 
on 22 May, 2009 with Magellan/IMACS as detailed in 
Green et al. (2010). 

For Chandra observations, we restricted our sam- 
ple to binary quasars with velocity differences Aw < 
800 km s^^. proper transverse separations Rp < 30kpc 
and redshifts z < 1. As we do not concern ourselves 
with binary quasars that do not meet these criteria in 
this paper, any such candidates can be considered as 
"discarded" for the piirposes of this work, although they 
will be published in a later paper (Myers et al., in prepa- 
ration). Our separation criterion selects hosts likely to 
be interacting on their first or second pass. The velocity 
criterion removes most chance projections but still allows 
for hosts in a variety of environments from isolated pairs 
to massive clusters. The redshift criterion prevents the 
necessary exposure times from becoming excessive. The 
properties of our final sample are shown in Tables 1 and 
2. 

3. Chandra X-RAY OBSERVATIONS 

We obtained X-ray images of the seven quasar pairs 
with the Chandra X-ray Observatory on the dates shown 

® uvxts=l OR lowzts=l 

® http://www.ucolick.org/~xavier/LowRedux/ 



in Table 1. We placed targets near the ACIS-S aimpoint, 
and tuned our exposure times to achieve ^100 counts for 
the fainter member of each pair, by converting the SDSS 
r mag to an expected fx using the 75th percentile X- 
ray-faintest value of log {fx / fr)=—0-5 from the Chandra 
Multiwavelength Project (ChaMP) QSO sample (Green 
et al. 2009). For every pair, we convert fx to ACIS-S 
counts/sec using PIMMS, with F = 1.9 through Ng"', 
and derive the exposure, which yielded exposure times 
from 12 to 30ksec, with a total of ~172ksec. 

3.1. X-rays from the Quasars 

The small (2 -3") separation of 3 of these pairs is not 
a challenge for Chandra. In all cases but one, the X- 
ray components are detected, well-resolved by Chandra, 
and correspond closely (< 0.2") to their SDSS coun- 
terparts. SDSS J160602.81-F290048.7 was not detected 
using wavdetect and a detection significance threshold 
corresponding to about one false source per ACIS chip. 
However, aperture photometry at the optical source po- 
sition shows 6 net counts, all above 2 keV. To avoid cross- 
contamination between QSOs in each pair, we extracted 
the X-ray photons from apertures corresponding to 90% 
of the counts (for 1.5 keV). Some of the QSOs in our 
sample yielded relatively few net counts. In such cases, 
instrumental hardness (photon count) ratios are often 
used, in the belief that genuine spectral fitting is not 
warranted by the data quality. There are several prob- 
lems with the use of hardness ratios (HRs). HRs do not 
take redshift or intervening Galactic hydrogen column 
into account. They are difficult to interpret or compare, 
since they convolve the intrinsic quasar SED with the 
telescope and instrument response, especially since the 
latter depends on both time and ACIS chip position^°. 
Finally, HRs waste useful spectral information by crudely 
binning counts. A direct spectral fit of the counts dis- 
tribution using the full instrument calibration, known 
redshift, and Galactic column Nh provides a much more 
direct measurement of qiiasar properties, most useful for 
comparison to other quasar samples. Correct spectral 
fitting does not underestimate errors; even in the low- 
count regime, one can obtain robust estimates of fit pa- 
rameter imcertainties using the Cash (1979) fit statistic. 
The spectral fitting we employ provides the most consis- 
tent and robust estimates of the physical parameters of 
interest — the power law slope and intrinsic absorption. 

We fit an X-ray power-law spectral model 

N{E) = N{Eo) exp[-7V^'^'a(£;)-7V^a(i?(l+z))] 

to the counts for each QSO using the CIAO tool Sherpa, 
where N{Eq) is the normalization in photons cm~^ sec~^ 
keV"-*^ at a reference energy £^0 (of 1 keV here), and a{E) 
is the absorption cross-section (Morrison & McCammon 
1983; Wilms et al. 2000). We fix N§''' at the appropriate 
Galactic neutral absorbing column taken from Dickey & 
Lockman (1990), and perform (1) a simple power-law fit 
with no intrinsic absorption component (model PL) and 
(2) include an intrinsic absorber with neutral column 
at the source redshift (model PLAbs)- Unbinned spectra 

See information on ACIS contamination at 
http: / /cxc. harvard.edu/cal 
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were fit using Cash statistics (Cash 1979). The best- 
fit model parameters for all components are shown in 
Table 2. 

The power-law energy index values T we measure are 
typical of SDSS Type 1 (broad line) QSOs in general 
(Green et al. 2009), with a mean of 2.14 ± 0.29 and 
a median of 2.11. Unabsorbcd fluxes and luminosities 
are calculated as detailed in the Appendix, using the T 
values from the PL^hs fits in every case except for the 
faintest object SDSS J160602.81+290048.7, where we as- 
sume r = 2.1. These values only differ substantially 
from the PL-only fit values in the two cases where there 
is absorption detected at >68% confidence. 

3.2. Comparison to Single Quasars 

Our small binary QSO sample has mean (median) red- 
shift 0.72 ± 0.18 (0.77). An excellent control sample is 
available already throiigh the ChaMP: wc have matched 
1175 SDSS QSOs from the SDSS photometric quasar cat- 
alog (Richards et al. 2009) to Chandra serendipitous X- 
ray sources measured in 323 X-ray images from Cycles 
1-6. To form a fair, high-quality comparison sample, we 
limit the ChaMP QSOs to those at z < 1.2, with expo- 
sure times > 4ksec, and off-axis angles 9 < 12'. This 
yields a control sample of 264 isolated QSOs, with mean 
(median) redshift 0.74 ± 0.32 (0.79), and a cumulative 
Chandra exposure of ~ 7.7 Msec. T?he ChaMP sample 
includes 70 QSO candidates (26%) with only photomet- 
ric redshifts (Weinstcin et al. 2004). Since quasars are 
known to have significantly higher fx //opt values com- 
pared to galaxies (e.g., Green et al. 2004), the additional 
criterion of X-ray detection for these photometric QSO 
candidates means that '^98% of them are indeed QSOs, 
as found in Green et al. (2009). Note that the par- 
ent SDSS QSOs were selected optically using the same 
approach as we used to target binary quasars. How- 
ever, since we specifically targeted the binary QSOs with 
Chandra, their X-ray data is of somewhat higher qual- 
ity (all on-axis, with a mean/median of 573/319 X-ray 
counts) compared to the control sample (135/42 counts 
in the mean). We find no significant difference in any 
of the measured ensemble properties. Comparing power- 
law fits, the mean (median) F is 2.14±0.30 (2.11) for the 
binary QSOs, and 1.96 ± 0.61 (2.01) for the comparison 
sample. 

To compare X-ray/optical luminosity ratios, as in 
Green et al. (2009), we first estimate the monochro- 
matic luminosity at 2500Aby finding the filter for each 
QSO whose de-redshifted effective wavelength centroid 
(taken from Fukugita et al. (1996)) is closest to 2500A 
in the restframe. We then asstune q;=0.5 for the 
optical continuum powerlaw slope to derive the rest- 
frame, monochromatic optical luminosity at 2500A in 
erg s~"'^Hz~^. Wc adopt the most common X-ray/optical 
measure for quasars, aox, defined as the slope of a hy- 
pothetical power-law from 2500 A to 2 keV i.e., aox = 
0.3838 \og{l^^^^ ^/l2kev)- The mean (median) aox is 

1.60 ± 0.21 (1.59) for the binary QSOs, and 1.57 ± 0.16 
(1.57) for the comparison sample. 

X-rays in quasars become weaker relative to optical 
emission as luminosity increases (Avni & Tananbaum 
1982; Wilkes et al. 1994; Green et al. 1995; Steffen et 
al. 2006; Green et al. 2009; Lusso et al. 2010). The 



Q^ox ( ^2500 A'' correlation is persistent across quasar sam- 
ples, but has large dispersion. One possible explanation 
of the observed trend is that AGN accretion may transi- 
tion between accretion states similar to those of Galactic 
X-ray binaries (XRBs), where different accretion rates 
harden or soften the overall SED. While AGN vary on 
much longer time scales than do XRBs, samples of ACN 
may show analogous trends in SED with luminosity to 
XRBs (Sobolewska et al. 2011). Statistical tests have 
shown that the ttox correlation is weaker with redshift, 
so that the aox ( ^2500 a-' '^^l^tionship is not a secondary 
effect of quasar evolution combined with the strong L — z 
trends of fiux-limited quasar samples. However, some 
dispute remains about the influence of selection effects 
(Yuan et al. 1998; Tang & Grindlay 2009). 

If our binary QSO systems genuinely reflect pairs at an 
unusual merging stage, perhaps being ignited or exacer- 
bated by an ongoing merger, we might expect to see dif- 
ferences in the properties of the AGN involved compared 
to a random selection of isolated quasars. In particular, 
if one or both nuclei in our pairs is being particularly af- 
fected by the merger, we might expect differences in the 
expected values of aox, given ^2-00 a ^'^^ each component 
of the pair. Figure 1 shows aox vs. optical 2500A log lu- 
minosity for the binary QSOs (black squares), with pair 
members linked by black lines. The comparison sam- 
ple of 264 z < 1.2 SDSS QSOs with Chandra detections 
from (Green et al. 2009) is also plotted, for which red tri- 
angles indicate spectroscopic redshifts, and blue circles 
show radio-loud objects. Binary QSO constituents ap- 
pear to follow the rather noisy trend of aox with optical 
luminosity. Only one QSO, SDSS J160602.81-F290048.7 
falls well away from the aox(? ° ) trend, with aox = 2.2 

2500 A 

at logZ ° = 30.39. This QSO is unusually faint in 

2500 A 

the X-ray band, and so may be a low redshift broad ab- 
sorption line quasar (BALQSO). This argument is fur- 
thered by the fact that it is probably X-ray absorbed (all 
6 measured photons are above 6keV). The vast majority 
of recognized BALQSOs in the SDSS are above redshift 
1.6 because only then is the CIV absorption rcdwards 
of the blue cut-off for SDSS spectroscopy.^^ In most 
cases, BALQSOs are X-ray weak due to large warm (ion- 
ized) absorbing columns (Green et al. 2001b; Gallagher 
et al. 2006). BALQSOs tend to have narrow H/3 broad 
line components, weak [OIII] lines, strong optical Fell 
emission — all of which are apparent in this object's SDSS 
spectrum — and be radio quiet. 
Criticism of the ensemble trend in aox( I :) ob- 

™^ 2500A; 

served in samples of isolated QSOs have been published 
(Yuan et al. 1998; Tang & Grindlay 2009), charging that 
it could be a selection effect caused by the different dis- 
persions in X-ray and optical luminosity, combined with 
flux-limited survey cutoffs. Binary QSOs represent ob- 
jects caught at the same epoch and in the same large- 
scale environment. If a sufficient sample of binary QSOs 
with dedicated followup also showed aox( I °) trends, 

2500 .A 

then this might obviate such criticism. However, since 
we are testing for potential effects of interaction be- 
tween the constituent QSOs, such a test is invalid. In- 

A much smaller number of the rare low-ionizatioii BALQSOs 
(with BALs just blueward of Mgll) are found at lower redshifts. 
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deed, we might test whether the differences between aox 
and ^ for each component of each pair (A aox and 
^^2500 A' '"^^P'^'^tivcly) are discrepant with the expected 
trends for isolated QSOs. As neitlier component of our 
pair is known to be special, we adopt a one-tailed dis- 
tribution and only allow these differences to be positive 
in value. We establish the background expectation for 
the relationship between A a^-^ and Al ° by select- 

2500 A 

ing 5000 pairs at random from the 264 SDSS quasars 
for which we have X-ray data from the ChaMP (note 
that there are then only N {N - l)/2= 34716 possible 
unique pairs to sample, so our precision cannot be in- 
creased without severely oversampling) . In Figure 2 we 
plot the distribution in density of our 35000 mock pairs 
in the (Auox, ^^2500 a-' P^^^® compared to the 7 data 
points. 

It is clear from Figure 2 that most of the pairs are not 
unusual as compared to background expectation. One of 
the data pairs is near the extreme of the distribution, 
with only 7% of mock pairs having similarly extreme 
values of (Aanx, Al °). However, as we are con- 

^ ° ' 2500A'^ ' 

sidering 7 pairs, a result at the 7% probability level is 
not unusual — indeed it should be expected. We demon- 
strate this further in the right-hand panel of Figure 2, 
for which we draw 5000 sets of 7 pairs at random from 
our 34716 unique mock pairs and plot the contour (from 
Figure 2, i.e. the 7% quoted in this paragraph) of the 
least likely pair in the distributions of 7 mock pairs. The 
histogram in Figure 2 demonstrates that most random 
sets of 7 pairs have one pair with a probability at the 7% 
level. Our data are therefore not unusual in the (Actox, 
^^2500 A-' P^^"^^' suggesting that either these values are 
not unusual for activated nuclei in ongoing mergers, or 
that we are not seeing a set of 7 ongoing mergers on the 
data. 

4. NEAR INFRARED IMAGING 

To search for extended host galaxy emission and/or 
morphological signs of mergers or interactions, we pro- 
posed near-IR imaging to optimize the contrast between 
the relatively blue quasar point source emission and the 
stellar light from the host galaxies. We were awarded 
2 nights to image binary QSOs on Mt Hopkins using 
the SAO Wide-field InfraRed Camera (SWIRC) on the 
6.5m MMT. SWIRC has 2048 x 2048 pixels spanning a 
^ 5.12' field of view with 0.15" pixels. We observed 9 
pairs from the larger binary sample in J (1.2 /im) band 
on the nights of 25 and 26 March 2010 and obtained be- 
tween 6 X 90 sec and 33 x 90 sec dithered images. We 
used the SWIRC pipeline to scale and subtract dark im- 
ages and remove sky background from all the images. 
The sky image per object frame was created using SEx- 
tractor (Bertin & Arnouts 1996). In the 5.12' x 5.12' 
field-of-view, each object frame contained at least three 
stars from the Two Micron All Sky Survey (2MASS) 
point source catalog (Cutri et al. 2003; Skrutskie et al. 
2006), which wc used to calibrate the astrometry of each 
frame and to determine the flux zeropoint from the mag- 
nitude conversions of Rudnick et al. (2001). The magni- 
tude where the number counts histogram turns over is a 
good general indicator of the limiting magnitude where 
incompleteness sets in. For our shallower field, a typi- 
cal exposure of 540 sec results in a limiting magnitude 



of 18.7 while for our deeper fields with exposure times 
of 2970 sec, the turnover magnitude is 20.3. We then 
used the imwcs software in the WCSTools package (Mink 
1997) to derive sky coordinates. We examined the distri- 
bution of FWHM for all images contributing to a given 
QSO field and excluded any outliers. We then stacked 
all the images of a QSO field using the Image Reduction 
and Analysis Facility (IRAF)^^ imcombine task — we av- 
eraged stacked science frames of all astrometrically cor- 
rected, sky-subtracted images, applying a Icr-clip. Small 
portions of the SWIRC field of view, especially the edges 
of each image, were disregarded due to significant con- 
tamination from CCD artifacts. Each of our stacked im- 
ages contains between 80 and 120 objects consistent with 
previous J-band surveys of the same depth (Saracco et 
al. 2001; Ryan et al. 2008). Seeing at the wavefront sen- 
sor varied from 0.7 to 1.1", yielding a typical median 
FWHM of 1.25"on our stacked images. 

We obtained SWIRC photometry for a total of four out 
of the seven Chandra-observed pairs in our sample. Near- 
infrared properties are given in Table 3. None of our 
SWIRC QSOs has 2MASS J-band counterparts, but we 
find excellent agreement between SWIRC and UKIDSS 
for the four QSOs with public UKIDSS photometry. We 
detect point sources for all QSOs, but no evidence of 
extended emission. Even SDSS J 1254-1-0846, a merger 
with spectacular tidal tails detected in our deep optical 
imaging (Green et al. 2010), shows no SWIRC evidence 
for disturbed morphology. 

5. QUASAR SPECTRAL ENERGY DISTRIBUTIONS 

To characterize SEDs, estimate bolometric luminosi- 
ties and check for starburst activity wc fit template SEDs 
to all QSO pairs in our sample. We fit to all of the near- 
infrared, optical and X-ray fluxes described above, using 
a library of 12 templates: a radio-quiet Type 1 quasar, 
a luminosity-dependent radio-quiet Type 1 quasar (i.e., 
where a„^ follows the known trend with I two 

°^ 2500 A^' 

Type 2 (narrow line) Seyferts, four starburst and four 
composite templates (Ruiz et al. 2010). We fitted all 
SEDs by using the minimization technique of Ruiz et 
al. (2007, 2010). Existing optical spectroscopy for each of 
our sources removes the uncertainty of using photometric 
redshifts for the SED fitting and provides a direct testbed 
for the accuracy of the fitting method. Broad-band SEDs 
for all 14 sources are given in Figure 3. Table 4 gives the 
different parameters of our best SED fits. 

In all 14 cases, cither a radio quiet Type 1 QSO 
(Richards et al. 2006 for u > 10^^ jjz and' Elvis et al. 
1994 for v < 10^^ Hz) or an AGN luminosity depen- 
dent template (Hopkins et al. 2007a) is needed to fit 
the photometry, consistent with the fact that each of 
our sources are spectroscopically identified as broad line 
QSOs. Nonetheless, for at least one component in 5 out 
of 7 pairs — 8 out of 14 sources in total — the best-fit SED 
requires an additional starburst component. 

We can contrast this fraction (8/14, or 57%) with that 
for single broad-line AGN (BLAGN) with spectroscopy 
from the ChaMP (Green et al. 2004; Green et al. 2009). 

IRAF is distributed by the National Optieal Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under the cooperative agreement with 
the National Science Foundation. 
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Of 758 spectroscopically identified ChaMP BLAGN, f 84 
have near-IR photometry, and 98 of those (54%) require 
a template with a starburst component. If we further 
restrict the ChaMP sample to the same 0.44 < z < 1 
redshift range as our binary QSOs, the fraction does not 
change (53%). To summarize, our binary QSO SEDs 
are no different than those of X-ray-selectcd QSOs with 
optical and near-IR counterparts in the ChaMP. 

When the AGN and/or starburst component contri- 
bution is estimated over the 10^^ - 10^® Hz wavelength 
range, where we have available photometry, the luminos- 
ity of all 14 sources appears to originate mainly from 
an AGN component (>55%). However, in the case of 
the 8 sources that require a starburst component, star- 
formation activity contributes at least 20% of the lumi- 
nosity emitted between 10^^ - 10^^ Hz. When we in- 
tegrate luminosity over the entire radio to X-ray wave- 
length range, starburst activity becomes the dominant 
component (>90%) in 6 cases, which may be indicative 
of intense star-formation events in their hosts. We warn, 
however, that the bulk of the starburst template contri- 
bution comes from longer wavelengths — far-IR to radio — 
than our available data. 

To seek further constraints on the mid- to far-infrared 
spectral regions, we have cross-correlated our binary 
QSO sample with the Preliminary Data Release cata- 
log (April 14, 2011) of the Wide-field Infrared Survey 
Explorer (WISE; Wright ct al. 2010). The catalog covers 
about 23,600 deg^ with typical stacked exposures near 
100 sec for at 3.4, 4.6, 12, and 22 /xm with an angular 
resolution of 6.1", 6.4", 6.5" and 12.0", respectively. Al- 
though WISE is thus unable to resolve the two compo- 
nents of our pairs, our systems are bright enough for 
WISE to be able to detect the emission from the paired 
system. Four of our pairs (J0740-h2926, J0813-F5416, 
J1508-F3328, J1606-F2900) have WISE counterparts with 
detections in all four WISE bands. To test our SED 
fitting results at these longer wavelengths, we have uti- 
lized the WISE detections to create the combined "pair" 
SEDs of the four pairs with mid-infrared detections. For 
each pair, we have simply added the near-infrared, op- 
tical and X-ray fluxes together for the two binary QSO 
components and appended the WISE fluxes into the to- 
tal pair SED. We then fit the summed SEDs for each 
pair. In the case of J0740-I-2926 and J1508-t-3328, the 
best solution is a pure AGN. In the case of J0813-I-5416 
and J1606+2900, the best solution is a composite SED 
with the AGN contributing up to 80% and 87% of the 
emission, respectively with the remainder coming from a 
starburst component. The agreement between these re- 
sults from the individual and paired fitting is almost per- 
fect. Disagreement for the case of J1508+3328, with no 
starburst component required in the summed SED with 
WISE, is likely due to the fact that it has the smallest 
total starburst contribution of all summed pairs (<12%; 
see Table 4). 

The predicted SED fits suggest that several of the 
QSO pairs in our sample may have significant ongoing 
starformation, detectable even in the presence of lumi- 
nous QSO emission. Interestingly, the one system, SDSS 
J1254-F0846 (Green et al. 2010), that is known to in- 
habit a merger, does not require a significant contri- 
bution from star formation. Previous studies of X-ray- 



selected (Trichas et al. 2009) and spectroscopically con- 
firmed (Lutz et al. 2008; Trichas et al. 2010; Kalfount- 
zou et al. 2010; Symeonidis et al. 2010) QSOs with far- 
infrared detections have shown that the vast majority 
of these sources are composite objects with very strong 
ongoing starburst events. While these studies were fo- 
cused on the brightest and rarest examples, siibsequent 
studies of submillimeter-detected Type 1 QSOs (Lutz et 
al. 2010; Hatziminaoglou et al. 2010) have made it clear 
that the sub-millimeter colors of Type 1 QSOs are sim- 
ilar to those of star-forming galaxies. This hints at an 
emerging picture where star formation is present in the 
environs of all AGN, consistent with merger models like 
those of Hopkins et al. (2005). On the other hand, in 
the local Universe, all black hole accretion as detected 
by hard X-rays is strongly disassociated with star forma- 
tion implying that there is a fundamental anti-correlation 
between the two that is not a selection effect (Schawinski 
et al. 2009). In the latter case, the prediction of starburst 
activity in the majority of the QSO pairs in our sample 
has strong implications for the dynamics of these merger 
systems that should be further investigated. 

6. SEARCH FOR EVIDENCE OF HOST CLUSTERS 
6.1. X-rays from Host Clusters 

Despite the presence of bright quasar emission, we 
know that Chandra can detect extended cluster emission 
against a typical ACIS-S background with '~50 diffuse 
counts or more in any of our fields (Green et al. 2002; 
Aldcroft & Green 2003). To ensure that we could de- 
tect clusters as weak as ^ 0.1 L'^, we slightly increased 
our proposed Chandra exposure times above what was 
required for the QSOs themselves (see §3 where nec- 
essary, extrapolating from typical cluster relationships 
(io.5-2feey^3 X lO''-''; Mullis et al. 2004). For two pairs 
we thus increased exposure times slightly: SDSSJ0740 
(-F5ksec) and SDSSJ1606 (+6ksec). 

The ACIS image of SDSS J1158-I-1235 (obsid 10314) 
displays significant extended X-ray emission — but, it is 
43" SSW of the QSO pair's midpoint. The peak of 
the extended X-ray emission is coincident with a lumi- 
nous i = 17.18 absorption-line galaxy at z = 0.2652 
(SDSS J115821.96-I-123438.6). At absolute magnitude 
Mi ^ —23.71, this is clearly the cD galaxy of an X- 
ray cluster. Using a circular aperture of 24" radius for 
the cluster, and a background annulus from 62 - 110" 
excluding all detected source regions, we derive 301 ± 19 
counts from the cluster. Assuming a Raymond-Smith 
plasma with T = 2keV and metallicity 0.2 solar, we 
derive using the Chandra Portable, Interactive Multi- 
Mission Simulator PIMMS^^ /(0.5-2keV) = 3.97 x 
10^^'* erg cm~^ s~^, and a luminosity of 8.2x 10**^ erg s""'^. 
Errors on these values are dominated by the spectral as- 
sumptions, and are likely to be ~15%. 

Otherwise, no significant extended emission sources are 
evident to the eye on the ACIS-S3 images in the imme- 
diate neighborhood of the QSO pairs. When searching 
for faint extended sources, however, it is important to 
minimize background contamination. The ACIS parti- 
cle background increases significantly below 0.5 keV and 
again at high energies. To optimize detection and visual 

http://cxc.harvard.edu/toolkit/pimms.jsp, originally Mukai 
(1993). 
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inspection of possible weak cluster emission, we first fil- 
tered the cleaned image to include only photons between 
0.5 and 2 keV. Around the QSO positions as detected by 
wavdetect, we masked out pixels within twice the ra- 
dius that encompasses 95% of the encircled energy. (The 
95% PSF radius at 1.5 keV is about 2.06".) For visual 
inspection, we also excised regions around all other de- 
tected sources, corresponding to 4 times the 4ct Gaussian 
source region output of wavdetect. 

We then generated 7 annuli of 50 kpc projected width 
each, starting at i? = 75 kpc from the mean of the de- 
tected QSO coordinates.^^ Though the sample redshifts 
range from 0.44 to 0.978, these radii only differ slightly 
between the targets (dispersion in the mean is about 
13%), so we used a single set of six 7" annuli from 10 
to 52". We set a background annulus from 60 - 110", 
and calculated radial surface-brightness profiles. 

There are just 2 fields with radial profiles that 
rise consistently inward towards the QSOs. For 
SDSS J0740-I-2926 (obsid 10312), the profile arises from 
some faint diffuse emission with a centroid about 7.5" 
NW of the mean QSO positions. The emission only en- 
compasses about 9.8 ± 3.6 net (0.5-2 keV) counts, and 
there is at least one other such source in the field, so 
we discount its reality. The other field with a suggestive 
radial profile is that of SDSS J1508-h3328 (obsid 10317), 
which similarly shows an apparent weak diffuse emission 
region at 6.8" W of the mean QSO position, with 14 ± 4 
net (0.5-2 keV) counts. There is no other such source ap- 
parent in the field. These weak excesses may represent 
the emission from nearby galaxies that fall individually 
below the detection level, or from a weak ICM 

The ACIS image of SDSS J0813+5416 (obsid 10313) 
shows signs in the smoothed image of extended emis- 
sion that could be more filamentary in shape, and so 
would not register as a significant trend in a radial pro- 
file plot. The emission appears to extend about 80" from 
SE to NW. Excluding the QSO regions, and using el- 
liptical source and background apertures (of about 0.7 
sq. arcmin area), we tally 64 ± 12 net source counts. 
There are no obvious optical counterparts that might 
be galaxies associated with a cluster merger or cosmo- 
logical filament. Assuming a Raymond-Smith plasma 
with T = 2keV and metallicity 0.2 solar, we derive 
/(0.5-2keV) = 7.43 x lO'^^erg cm'^ s'^. If at the 
z = 0.779 redshift of the QSOs, the cluster luminosity is 
2.0 X 10^^ erg s'^ 

6.2. Optical Imaging 
6.2.1. Kitt Peak 

To study each component of each binary quasar in the 
optical, and to search for further signs of merger activ- 
ity or local galaxy overdensities, we imaged six binary 
QSOs at Kitt Peak National Observatory (KPNO) us- 
ing the 4m Mayall telescope on the nights of 2009 March 
17-19. All images were acquired with the MOSAIC 8K 
camera (8192 x 8192 pixels; 0.26" pixel"^) in one or more 
filters using the r, i, and z bandpasses. Integration times 
ranged from 900 to 9000 s, depending on the filter and the 
redshift of the target binary quasar. The seeing varied 

1"* Since SDSS J1606-I-2900B was not detected by Chandra, we 
simply use its optical position from SDSS. 



during the observing run from 0.77" to 1.69" (FWHM), 
as measured from the combined frames. 

Image reduction was conducted using the mscred pack- 
age within IRAF. Processing of the raw images involved 
the standard procedure of bias correction and flat-fielding 
using dome flats and deep sky exposures. After initial 
processing, individual images were astrometrically cor- 
rected and median combined to yield a higher S/N im- 
age. 

Object detection and photometry was conducted us- 
ing SExtractor (Bertin & Arnouts 1996) via the ChaMP 
image reduction pipeline (Green et al. 2004). Since all 
images were acquired during non-photometric sky con- 
ditions, instrumental magnitudes were transformed to 
the standard system by calibrating to overlapping SDSS 
DR7 data using dereddened magnitudes. We typi- 
cally achieved magnitude limits'^^ of 24 24.5 in i and 24 
in r. We examined all images for any evidence of ex- 
tended emission or disturbed morphology. This led to 
the discovery of tidal tails in both r and i band images 
of SDSS J1254-I-0846. However, due to poor weather, 
we were only able to image 6 of the 7 fields in this 
Chandra subsample (all but SDSS J1418+24410), and 
obtained imaging in more than one band for only 2 fields: 
around SDSS J1158-hl235 and SDSS J1254+0846. This 
precluded an effective photometric search for galaxy over- 
densities, described below. 

6.2.2. SDSS 

With optical imaging of adequate depth, we can pho- 
tometrically detect an overdensity of galaxies — because 
early-type galaxies at a given redshift have a narrow 
range of colors which form a cluster "red-sequence" 
(Gladders & Yee 2000) in their color-magnitude diagram 
(CMD). In the neighborhood of a QSO pair, the most 
convincing optical cluster detection would have a large 
number of optical galaxies clustering at small distances 
from the QSO pair center, and those galaxies would have 
well-measured colors clustering at small distances from a 
single locus in the CMD. We therefore define a distance- 
and error- weighted color mean (DWCM), given by 

DWCM = 5:l!:^/^^, (1) 

where cr^ = cr?^_,N -I- i??, and R is the projected distance 
from the mean quasar pair position in units of Mpc at 
the QSO redshift. Thus a bright galaxy with a small 
color error could contribute as much to the DWCM as 
a fainter galaxy closer to the center point. With the 
IMpc scafing, a 250 kpc projected galaxy distance and 
a typical color error of 0.25 contribute about equally to 
the weighting. Using the same DWCM calculation for 
any number of randomly-chosen locations in the same 
large-field optical image of our quasar field allows us to 
quantify the significance of the DWCM measured around 
our QSO pairs, in a way that naturally accounts for the 

15 We compare SDSS modelJlag to SExtractor MAGJVUTO val- 
ues. 

1*' Wc quote the magnitude where the number counts peak in a 
differential (0.25 mag bin) number counts histogram. This corre- 
sponds approximately to 90% completeness in the magnitude range 
20-25, and is typically about 1 mag brighter than the 5cr limiting 
magnitude (Green et al. 2004). 



8 



Green et al. 



characteristics of the relevant imaging such as depth and 
image FWHM. 

We selected all objects within 48' of each quasar pair's 
mean position from the SDSS DR8 using the CasJobs 
query interface. The turnover (model) magnitudes are 
r ~ 22.8, r ~ 22.2 and i ~ 21.4, which correspond 
to about 50% completeness for point sources^^. We in- 
cluded only objects with 15 < r < 22.5, for which the 
median error in (r — i) and {i — z) color are 0.105 and 
0.192, respectively We calculated the DWCM in both 
(r — i) and {i — z) for 1000 random positions for each of 
our quasar pairs, and at the actual position of the quasar 
pair, always in an annulus between 25kpc and 500 kpc 
projected radius at the QSO redshift. We then com- 
pare the DWCM to the expected observed-frame color 
of the galaxy red-sequence for the Schechter magnitude 
Ad* based on the redshift of each quasar pair, adopting 
the red sequence models of Kodama & Arimoto (1997) 
transformed to the SDSS filters (T. Kodama 2004, pri- 
vate communications). 

If the DWCM at the QSO position is appropriate for 
the red-sequence color expected at the pair redshift, and 
the red-sequence scatter is small, this could be especially 
convincing evidence for a physical galaxy cluster. To es- 
timate the prominence of the red-sequence, we calculate 
the variance in DWCM and compare the results for our 
set of random locations with results for the location of 
the quasar pair. 

Only the SDSS J1158-hl235 position shows a DWCM 
both close to the (i — z) color expected for an overdensity 
and significantly different from contours for a DWCM 
derived by sampling random positions in the field. How- 
ever, this test does not hold for r — i. Both plots are 
shown in Figure 4. 

We conclude that we detect no significant galaxy den- 
sity enhancements of the color and magnitudes expected 
for early-type galaxies at the redshift of our QSO pairs. 
However, we note that across the full redshift range of 
our sample, the expected apparent Schechter magnitude 
in the r band ranges from 20.8 at z = 0.44 to 23.9 at 
z ^ 1. Therefore, for objects beyond z ~ 0.7, we detect 
few bright galaxies of the expected color. 

7. DISCUSSION 

In binary quasars, both host galaxies would already 
have substantial SMBH and pre-existing stellar bulges, 
so must have a significant history of accretion before the 
observed episode of simultaneous activity. But at any 
moment in time (i.e. when observed) either or both com- 
ponent QSOs might otherwise be quiescent. Our intent 
in this work is to study both the SEDs and environ- 
ments of binary quasars to search for signs that inter- 
action might indeed be triggering the currently observed 
activity. One alternative to the interaction/triggering in- 
terpretation is simply that QSOs are more likely to be 
found in overdense regions, with a QSO pair likely to be 
found in some fraction of those, perhaps more likely in 
those inhabiting the largest overdensities. 

To probe for host signatures of merging or triggering is 
challenging in luminous QSO pairs, both because of their 
bright nuclei, and because they are found at significant 
redshifts, making host imaging difficult. Nevertheless, in 

Based on comparison to http://www.sdss.org/dr7/products/ge 



our small Chandra sample of 7 binary quasar pairs, we 
have discovered one clear example of an interacting sys- 
tem in our lowest-redshift pair, SDSS J1254-I-0846 (Green 
et al. 2010). In this paper we pursued two other potential 
indicators of unusual accretion or star-formation activ- 
ity: multiwavelength spectral energy distributions and 
environment. 

7.1. Spectral Energy Distributions 

Analyzing results from published optical/infrared pho- 
tometry and our own Chandra observations, we find that 
the SEDs are consistent with those of isolated QSOs. 
Their X-ray spectra are typical, and show no sign of ex- 
cess absorption that might be expected in systems with 
accretion rates enhanced by interactions that dissipate 
angular momentum of gas. The ratio of optical to X- 
ray emission in these QSOs, characterized here by ctox is 
also typical, both in its distribution and in its correlation 
with luminosity. Such a finding might be expected based 
on these pairs' original selection by their typical optical 
quasar colors Myers et al. (2008). 

Based on our SED fits, the available optical and 
near-IR SEDs show possible evidence for enhanced star- 
formation activity, because the best fit requires a star- 
burst template in addition to a standard QSO tem- 
plate. It would be of great interest to test more robustly 
whether this tendency is statistically different from iso- 
lated QSOs. In a subsequent paper (Trichas et al. 2011), 
we are planning to utilize the large number of ChaMP 
spectroscopically identified isolated QSOs to compare 
their SEDs to a much larger sample of spectroscopically 
identified QSO pairs (e.g., Myers et al. 2008). Inclusion 
of WISE and/or Spitzer, Herschel and ALMA photom- 
etry would greatly improve over current constraints on 
star- formation activity. 

For a larger binary quasar sample, we can also correlate 
the SED characteristics and LboI with dynamical charac- 
teristics like Rp and Av — do smaller separations and/or 
lower velocities result in more luminous, high column 
systems? 

Today's favored models (e.g., Hopkins et al. 2008) as- 
sociate luminous AGN activity with major mergers, so 
the lack of significant SED-based evidence for interac- 
tions is interesting. Proximity does not dictate merg- 
ing. Differences between predicted galaxy-galaxy merger 
rates can be significant (factor ~ 5; Hopkins et al. 2010), 
attributable at least in part to the treatment of the bary- 
onic physics, especially those in satellite galaxies. 

Over a redshift range similar to our sample, Cisternas 
et al. (2011) find no significant difference in the frac- 
tion of (HST ACS) distorted morphologies between X- 
ray active and inactive galaxies in the COSMOS field. 
Consequently, they argue that the bulk of black hole ac- 
cretion has not been triggered by major galaxy mergers, 
but more likely by alternative mechanisms such as inter- 
nal secular processes or minor interactions. 

7.2. Environments 

We find no evidence that these pairs inhabit significant 
galaxy overdensities based on a search for red-sequence 
galaxies in SDSS optical imaging. Neither do they show 
signs of inhabiting a hot ICM that might be associated 
with a significant cluster of galaxies or a massive dark 
teatt)Jrt4iftte.html 
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While wc might hope that quasars — especially binary 
quasars — would be signposts for high-rcdshift clusters, 
this has not turned out to be the case. At low redshift, 
there are just a handful of X-ray clusters associated with 
quasars or powerful radio galaxies at lower redshifts (e.g., 
Cygnus A: 3C 295, Allen et al. 2001; IRAS 09104+4109, 
Iwasawa et al. 2001; HS1821+643, Russell et al. 2010), 
and even fewer at high redshifts (Siemiginowska et al. 
2010). 

The fraction of galaxies hosting AGN evolves with cos- 
mic time (Shi et al. 2008; Martini et al. 2009; Shen et al. 

2010a; Haggard et al. 2010) and is likewise affected by 
environment (e.g.. Strand et al. 2008). Luminous quasars 
and intense star formation activity both tend to be found 
at z> 1, where there are still very few massive clusters 
known. While the space density of luminous AGN de- 
creases drastically toward the present day (e.g., Silver- 
man et al. 2008), the clusters are assembling. Burgeon- 
ing detections of galaxy clusters based on the Sunyaev- 
Zeldovich effect (SZE) (Staniszewski et al. 2009; Vander- 
linde et al. 2010) may help widen the overlap. 

Whether AGN favor or eschew cluster environments in 
a given epoch is another question. At low redshifts, the 
fraction of galaxies that host X-ray AGN appears to be 
the same in clusters as in the field (Haggard et al. 2010), 
although the fraction in clusters may evolve more rapidly 
than the field (Martini et al. 2009). 

For galaxies at low redshift [z < 0.1), lower density 
environments have fractionally more galaxy pairs with 
small projected separations and relative velocities (Elli- 
son et al. 2010). Conversely, selection of pairs by small 
projected separation and low Av tends to select lower 
density environments, an effect which may apply here as 
well, since we restricted the parent binary quasar sample 
to velocity differences Av < 800kms~^ and separations 
Rp < 30kpc. Galaxies in the lowest density environ- 



ments show the largest star formation rates and asym- 
metries for the smallest separations, suggesting that trig- 
gered star formation is seen only in lower density envi- 
ronments (Ellison et al. 2010). Whether this does or does 
not apply to AGN triggering remains unclear. 
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APPENDIX 

APPENDIX: EXPLICIT X-RAY FLUX AND LUMINOSITY CALCULATIONS 

We often assume the monochromatic flux density for an underlying intrinsic power-law to have form f oc v", where 
/ is the monochromatic flux (e.g., in erg cm~^ s~^Hz~^) and v is the power-law frequency index. For X-rays, the 
photon index F is more commonly used, where a = (1 — F). 

We flt an X-ray power-law spectral model 

NiE) = NiEo) (§-) eM-N^ME) - NI,a{E{l + z))] 

to the X-ray counts as a function of energy, where N{Eo) is the normalization in photons cm~^ sec~^ keV~^ at a 
chosen reference energy Eq, F is the photon index, and ct{E) is the absorption cross-section. We flx N§°'^ at the 
appropriate Galactic neutral absorbing column, and allow for an intrinsic absorber with neutral column at the 
source redshift. 

The X-ray monochromatic energy flux without the effects of absorption is 

/ E \(i-r) 
f{E)=EN{E) = EN{Eo) (— ) 

in keV cm~^ sec~^ keV~^. Then, since f{Eo) = Eq N{Eo), we can express the monochromatic energy flux as 

/ E \(i-r) 
f{E) = fiEo) (— ) . 

To obtain the more standard units of erg cm^^ sec^ Hz^^, multiply by 6.629 x 10"^'' (from conversion factors 1.602 x 
lO-^erg/keV and Bz-^=hke\-^ where h = 4.138 x lO'^^ keV/Hz). 
The integrated flux observed between energies Ei and E2 is 

F= r mdE= [Et^^ - e['-^^] _ N{Eo) [Et^^ - Ei'-^^] 



El 



E, 



(i-r) 



(2-r) 



E^ 



(2-r) 
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If above is in units of keV cm ^ sec ^, multiplying by 1.602x10 ^ yields observed broadband flux in erg cm ^ sec^. 

Note that as F — > 2, via '. 
(or luminosity) to another 



Note that as F — > 2, via L'Hospital's rule F — > ^^1^ ln{E2/E\). Note also that to convert from one broadband flux 



FjE^ -Ej) _ [E: 
F{E, - E^) [E, 



(2-r) 



E. 



(2-r)i 



(2-r) 



E. 



(2-r)i 



Due to the redshift, the measured spectral flux /j, is related to the spectral rest- frame luminosity L^i , where 
z/' = (1 + z) V, as 

_{l + z) L^, 

Jv — 



The factor of (1 + z) accounts for the fact that the flux and luminosity are not bolometric, but are densities per unit 
frequency. (The factor would appear in the denominator if the expression related flux and luminosity densities per 
unit wavelength.) 
The monochromatic luminosity is therefore 



4 7rd| _ 47rd| 



(1 + ^) 



but since fi, oc and a = (1 — F), 



so that 



(1 + ^)-"^ V/., 

(i-r) 



V =47rdi(l + z)(r-2)/^, 



in erg sec~^ Hz~^. In this way, the flux measured at v in the observed frame yields the monochromatic luminosity L^' 
in the rest frame. 
The broadband luminosity in erg sec~^ is therefore 



Lx= \ i(z/)rf!/ = 47rrf|(l + ^)(^-2) / f{y)dv. 



Then, since 



J 1/1 J 1/1 



(-) 



(i-r) 



/(^o) 
„(i-r) 



/2-r) 



(2-r) 



(2-r) _ (2-r) 

^2 ^1 



,(i-r) 



we get 



T — /I --^2 n I ~\(r-2) fi^o) 
Lx = ^T^di^{l + zy (Try 



(2-r) 



f2-r) „(2-r) 



(2-r) 



where the final term is convenient for L'Hospital's rule. Perhaps more intuitively, we can write 

L. = 4.4(l + ^)(-^^ [^^^^^^^^^^]. 
To substitute E in keV for frequencies above, just multiply by 2.41666 x 10^^ Hz/keV. 
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Fig. 1. — : aox vs. optical 2500A log luminosity for the binary QSOs (black squares), with pair members linked by 
black lines. The comparison sample of 264 C/iandra-detected SDSS QSOs with 2; < 1.2 from Green et al. (2009) is also 
plotted, for which red triangles indicate spectroscopic redshifts, and blue circles show radio-loud objects. 
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Fig. 2. — : LEFT: Distributions of A a^^, Al^^^^j^ for individual quasars and our 7 pairs of close quasars. The contours 

are the density of 5,000 pairs of quasars drawn at random from the 264 SDSS quasars in the ChaMP. Points are the 7 
genuine pairs of quasars discussed in this paper. The most extreme of the 7 data points has a 7% probability of being 
drawn at random from the distribution of possible pairs of quasars. RIGHT: As our sample represents drawing 7 pairs 
of quasars, rather than just one pair, we repeat our experiment but testing instead the most improbable pair drawn 
at random in 5,000 samples of 7 mock pairs. We histogram the value of the contour (i.e. from the left-hand panel) for 
the most improbable of the 7 mock pairs. The dashed line is for the actual binary quasar sample data. 
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Fig. 3. — : Near-infrared to X-ray SEDs (Ruiz et al. 2010) in vfi, for each of our quasar pairs. Solid blue lines show 
the total predicted SED. Green and red lines are the corresponding AGN and starburst templates used, respectively. 
Parameters for model fits are given in Table 4. 
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Fig. 3. — : Near-infrared to X-ray SEDs (Ruiz et al. 2010) in vfi, for each of our quasar pairs. Sohd blue hues show 
the total predicted SED. Green and red lines are the corresponding AGN and starburst templates used, respectively. 
Parameters for model fits are given in Table 4. 
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Fig. 3. — : Near-infrared to X-ray SEDs (Ruiz et al. 2010) in vfi, for each of our quasar pairs. Solid blue lines show 
the total predicted SED. Green and red lines are the corresponding AGN and starburst templates used, respectively. 
Parameters for model fits are given in Table 4. 




Variance Variance 

Fig. 4. — : Test of optical cluster significance. (LEFT) Contours from 10 to 90% in relative frequency show the distance- 
weighted (i — z) color measure (DWCM) vs. its variance for 1,000 random positions within the (2deg) SDSS field of 
the QSO pair SDSS J1158-I-1235, calculated using the DWCM prescription described in § 6.2.2. The asterisk marks the 
DWCM and its variance calculated at the actual position of the quasar pair SDSS J1158-I-1235. The expected location 
{i — z = 0.98) of the galaxy red sequence at to* based on the mean redshift z ^ 0.878 of the quasar pair is shown as 
a solid horizontal line, adopting the red sequence models of Kodama & Arimoto (1997). While the variance is typical 
80%) for the field, the actual mean color (DWCM) is both unusual and close to the expected red-sequence color. 
(RIGHT) The DWCM vs. its variance for (r — i) shows that the expected red-sequence value at the midpoint of the 
QSOs' positions (asterisk) is entirely typical for random positions in the field, and far from the expected red-sequence 
color (horizontal line). 



TABLE 1: 
Chandra Close Binary Quasar Sample 



Pair Name ObsID Exposure 6 ^ Rp^ Galactic Ntt'^ 

(sec) (arcscc) (kpc) (10'^" cm ) 



SDSS J0740+2926 


10312 


20859 


2.6 


15.0 


4.24 


SDSS J0813+5416 


10313 


30625 


5.0 


26.9 


4.21 


SDSS J1158+1235 


10314 


30827 


3.6 


17.0 


2.07 


SDSS J1254+0846 


10315 


15967 


3.8 


15.4 


1.92 


SDSS J1418+2441 


10316 


29762 


4.5 


21.0 


2.00 


SDSS J1508+3328 


10317 


31317 


2.9 


16.0 


1.51 


SDSS J1606+2900 


10318 


12852 


3.5 


18.4 


3.19 



^Separation between QSO components in arcsec. 
''Proper separation between QSO components in kpc. 

Galactic column in units 10'^'^ cm^'^ from the NRAO dataset of 
Dickey & Lockman (1990). 



TABLE 2: 
Chandra Close Binary Quasar Sample 



SDSS NAME 




r mag 






Counts'^ 




1 Pi 


i PLAbs 




log/x^ 


log/x* 




J074013.42+292645. 


.7 


19.47 


0, 


.978 


H 


82.4 


1, 


oK+0.20 

•°>^-0.20 


1 oc+0.26 
1.00_Q 20 


<27 


-13.527 


25.993 


1.583 






1 Q 07 


u. 


.9803 


S 


^00.0 


z. 


94+0.11 
•^4-0.11 


9 94+0.11 


<y 


1 Q r\A A 


zD.ooy 


1 ^An 


J081312.63+541649. 


.8 


20.08 


0, 


.7814 


K4m 


195.7 


2, 


CI +0.15 
•^^-0.14 


9 C1+0.15 
^•'->1_0.14 


<12 


-13.414 


26.006 


1.426 




Q 


1 7 1 S 


n 

u. 


.7792 


S 


ZOOD.O 


1 

1. 


qf,+0.06 
•="-'-0.06 


, qo+0.04 

1.9Ci_g Q4 




-IZ.ZOO 


97 IT^I 
Z i .UOl 


1 Af{C\ 


J115822.77+123518, 
J115822.98+123520. 


,5 
,3 


19.85 
20.12 


0, 
0, 


.5996 
.5957 


M08 
M08 


367.7 
413.6 


2, 
2, 


9f;+0.10 

ic+O.OQ 
■^"-0.09 


9 c;i +0.19 
^■''-'--o.ig 

^■-'-^-().14 


<9 


-13.098 
-13.069 


26.017 
25.999 


1.335 

1.292 


J125454.86+084652. 


,1 


19.43 


0, 


.4401 


G 


349.5 


2, 


1 1 +0.10 
-'•^-0.10 


9 1 -1 +0.10 
^•^^-0.10 


<2 


-12.853 


25.887 


1.355 




Q 


1 7 08 
1 / .uo 


n 

u. 


.4392 


G 


J. 1 yo.o 


2 


•"*-0.04 


9 (14+0.04 




-12 129 




1 431 


J141855.41+244108. 


,9 


19.21 


0, 


.5728 


S 


864.9 


1, 


□4+0.06 
^^^-0.06 


1 99+"" 


<7 


-12.692 


26.305 


1.282 


J141855.53+244104. 


,7 


20.13 


0, 


.5751 


M08 


33.0 


0, 


79+0.28 
•'■^-0.28 


1 fi+0.60 
^■"-0.54 


1581113 


-13.769 


25.132 


1.575 


J150842.19+332802, 
J150842.21+332805. 


,6 
,5 


17.80 
20.19 


0, 
0, 


,8773 
.878 


S 
H 


1040.9 
111.1 


2, 
2, 


06 

■-'-^-O.Ofi 
07+0.45 
■■^'-0.40 


^•"^-0.20 


<8 
<27 


-12.670 
-13.708 


26.807 
25.860 


1.514 
1.479 


J160602.81+290048. 


.7 


18.35 


0, 


.7701 S 


6.1 


0, 


f;n+0.71 
■"^-0.73 


2.1 


<0.01 


-14.491 


24.844 


2.129 


J160603.02+290050. 


,8 


18.25 


0, 


.7692 M08 


144.1 


2. 


9n+0.43 
■^"-0.42 


2 44+0.17 


<6 


-13.171 


26.222 


1.611 



'^SDSS dcrcddcnod PSF magnitude. 

bRedshift. K4m - February, 2008 KPNO/4m; M08 - Myers et al. (2008); H - Hennawi et al. (2006); G - Green 
et al. (2010); S - SDSS 
'^Net 0.5-8 keV counts. 

''Best-fit X-ray power-law photon index. Uncertainties are the 68% confidence limits. If no uncertainties are 
shown, then the value is frozen to enable fitting of N^. For J160602. 81-1-290048. 7, we freeze F simply to fit the 
overall normalization. 

°Bcst-fit intrinsic column for PLAbs model in units 10^*^ cm~^. Upper limits are at 68% confidence. 
'X-ray flux (0.5-8 kcV) in erg cm^^ s^^ calculated using the PLAbs model. 
sX-ray luminosity at 2keV in erg Hz~i. 
^Uox, the optical/UV to X-ray spectral index. 



TABLE 3: 

Chandra Binary Quasar Sample Near-Infrared Properties 



SDSS Name 


Exposure 




err 




errj(7 " 


Y * 


errY ^ 


H " 


errH ' 


K j 


errK 


J074013.42+292645.7 
























J081312. 63+541649.8 

TOSI "^l 1 C\-\-^A^ fiAfi Q 
J uoioio. lUn^O'i-LU^u-y 


2610 


17.869 
18 744 


0.034 

U.UU-L 


















J115822.77+123518.5 
J115822.98+123520.3 


900 
900 


17.174 
17.411 


0.013 
0.016 










17.085 
17.515 


0.050 
0.074 


16.239 
16.732 


0.036 
0.056 


J125454.86-I-084652.1 
J125455.09-I-084653.9 


1260 
1260 


18.098 
16.063 


0.061 
0.010 


17.925 
16.079 


0.045 
0.009 


18.424 
16.505 


0.034 
0.008 


17.066 
15.433 


0.022 
0.006 


16.129 
14.327 


0.029 
0.007 


J141855.41+244108.9 
J141855.53+244104.7 
























J150842.19-I-332802.6 
J150842.21+332805.5 


630 
630 


16.551 
18.473 


0.058 
0.031 


16.629 
18.463 


0.012 
0.063 















J160602. 81+290048.7 
J160603.02+290050.8 



17.398 
17.369 



0.023 
0.022 



=-Total MMT-SWIRC exposure time in seconds. 
'•SWIRC J-band magnitudes. 

Error in SWIRC ,J-band magnitudes. 
'^UKIDSS J-band magnitudes. 
''Error in UKIDSS J-band magnitudes. 
fUKIDSS Y-band magnitudes. 
sError in UKIDSS Y-band magnitudes. 
hUKIDSS H-band magnitudes. 
'Error in UKIDSS H-band magnitudes, 
j UKIDSS K-band magnitudes. 
■^Error in UKIDSS K-band magnitudes. 



TABLE 4: Quasar Spectral Energy Distribution Fit Results 



SDSS Name 






Jsb'= 


PiRx" 




707401 S 42-l-2Q2fi4'i 7 
J074013.44+292648.3 


46.362 
46.909 


QSO 
QSO 




100 /o 
100/0 


100 /O 
100/0 


T0S1 'il 2 fiS-l-'i41 fi4Q 8 
J081313. 10+541646.9 


48.159 
47.087 


oso 
QSO 


M82 
NGC7714 


71/29 
79/21 


79/21 


J115822.77+123518.5 
J115822.98+123520.3 


47.185 
45.770 


LDQSO 
LDQSO 


IRAS 12112+0305 
NGC7714 


74/26 
58/42 


<10/>90 
58/42 


J125454.86+084652.1 
J125455.09+084653.9 


45.787 
46.457 


LDQSO 
QSO 




100/0 
100/0 


100/0 
100/0 


J141855.41+244108.9 
J141855.53+244104.7 


48.009 
47.713 


LDQSO 
QSO 


IRAS 12112+0305 
IRAS 12112+0305 


80/20 
76/24 


<10/>90 
<10/>90 


J150842. 19+332802.6 
J150842.21+332805.5 


47.013 
48.127 


QSO 
LDQSO 


IRAS 12112+0305 


100/0 
77/23 


100/0 
<10/>90 


J160602.81+290048.7 
J160603.02+290050.8 


48.225 
46.641 


LDQSO 
LDQSO 


M82 


68/32 
100/0 


<10/>90 
100/0 



''Log of the luminosity from 10^' - 10^®' Hz in units of erg s~^' from template fit. 

I'AGN template used in the best fit solution: QSO (radio quiet QSO), LDQSO (luminosity-dependent QSO template). 
"^Starburst template used in the best fit. 

"^Percent QSO/starburst contribution in the range 10^^' - 10^®' Hz. 
"Percent QSO/starburst contribution in the range 10^- - lO^''- Hz. 



